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ABSTRACT 



Context. The charged component of the insterstellar medium consists of atomic and molecular ions, electrons, and charged dust 
grains, coupled to the local Galactic magnetic field. Collisions between neutral particles (mostly atomic or molecular hydrogen) and 
charged species, and between the charged species themselves, affect the magnetohydrodynamical behaviour of the medium and the 
dissipation of electric currents. 

Aims. The friction force due to elastic collisions between particles of different species in the multi-component interstellar plasma is a 
nonlinear function of the temperature of each species and the Mach number of the relative drift velocity. The aim of this paper is to 
provide an accurate and, as far as possible, complete set of momentum transfer rate coefficients for magnetohydrodynamical studies 
of the interstellar medium. 

Methods. Momentum transfer rates are derived from available experimental data and theoretical calculations of cross sections within 
the classic approach developed by Boltzmann and Langevin for a wide range of values of the temperature and the drift velocity. 
Results. Accurate numerical values for momentum transfer rates are obtained and fitted to simple analytical formulae expressing 
the dependence of the results on the gas temperature and the relative drift velocity. The often used polarization approximation is in 
satisfactory agreement with our results only for collisions between H2 and molecular ions (HCO^, Hj). For other kinds of collisions, 
the polarization approximation fails by large factors, and must be replaced by more accurate expressions. 

Key words. Atomic processes ;Molecular processes;Plasmas;MHD;ISM:clouds;ISM:jets and outflows;ISM:magnetic fields 



1. Introduction 

The interstellar medium (ISM) is a multi-component plasma 
that consists mostly of hydrogen, helium, heavy ions, electrons 
and charged dust grains. The interaction between these compo- 
nents, and their coupling with the Galactic magnetic field, de- 
termine the dynamical properties of the ISM, control its evo- 
lution and the nature of the star formation process. In partic- 
ular, the momentum exchange in collisions between neutral and 
charged particles is responsible for transfering the effects of elec- 
tric and magnetic forces to the neutral component, allowing the 
magnetic field to drift out of weakly-ionized molecular clouds 
(Mestel & Spitzer 1956), damping the propagation of Alfven 
waves (Zweibel & Josafatsson 1983), and heating the gas by the 
frictional dissipation of turbulent energy (Scalo 1977). 

In a companion paper (Pinto, Galli & Bacciotti 2007, here- 
after Paper I), we have derived the equations governing the dy- 
namics of a three-fluid system, reducing the set of equations to 
a momentum equation for the mean fluid and an evolution equa- 
tion for the magnetic field, plus two relations for the drift veloc- 
ities in terms of the mean fluid velocity and the magnetic field. 

In this paper, we report on a detailed analysis of collisional 
rate coefficients involving the most abundant neutral and charged 
species in the ISM. The paper is organized as follows: in Sect.|2] 
we give the general expression for the friction force and the mo- 
mentum transfer rate coefficient for elastic collisions, and we ob- 
tain an analytical solution for a cross section varying as a power 



of the relative velocity; in Sect. [3] we consider collisions with 
H2 of HCO^, H3 , H^, and electrons, using available theoretical 
and/or experimental determination of the collision cross section; 
similarly, in Sect. |4] we consider collisions with H of C^, H^, 
and electrons; in Sect.|5] we consider collisions of H^ and elec- 
trons with He; in Sects. |6]and|7]we consider collisions between 
charged dust grains and neutral particles, and between charged 
particles, respectively; in Sect. [8] we give analytical approxima- 
tions for our numerical results; finally, in Sect.|9] we summarize 
our conclusions. 

2. Friction force and rate coefficients 

The general expression of the momentum acquired per unit time 
and unit volume ("friction force") by a particle of species s with 
mass Ills and initial velocity v, ("test particle") through collisions 
with particles of species s' with mass mj- and initial velocity Vj' 
("field particles") was given by Boltzmann (1896), 
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where /(v,) and f{\s') are the velocity distribution functions of 
the two species, Wj is the velocity of the test particle after the 
collision, dcr/dQ. is the differential scattering cross section, and 
Vss' is the relative velocity of the particles (before the collision). 
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For elastic collisions, the last term in Eq. ([TJ, representing 
the momentum change of the test particle after the collision, can 
be written as 
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(3) 



where fiss' - nisins' /(nis + nis') is the reduced mass of the system, 
dcr 



o"mt = I ^(1 -cos0)<iQ 



(4) 



is the momentum transfer cross section, and is scattering an- 
gle in the center-of-mass system. Eq. ([T]i reduces then to the ex- 
pression 

F,j' = yL/,j'n,n,' I d\J(\s) I £/Vj'/(v,Ov,j'Crn,t(v,' - v,). (5) 

The six integrations of eq. (|5]l require some care because of 
the presence of the relative velocity in the integrand. When the 
velocity distribution of both species is maxwellian with temper- 
atures r, and Ts', the integrations can be carried out explicitly, 
as shown first by Langevin (1905), and the result is 



F.„' = ass'iu,' -u.,), 



(6) 



where u.s and Uj' are the mean (maxwellian) velocities of each 
species and the friction coefficient ass' is defined as 



ass' = l^ss'nsns'{o-v}s 



(7) 



In this expression, {crv}ss' is an average of the momentum trans- 
fer cross section over the relative velocity of the interacting par- 
ticles, 

{o-v}ss' = ass' 
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where the nondimensional variables £, and x are defined by 
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and 

Vrf = \vis-W-s'\ 

is the drift velocity, and 



ass' = 



ZKl ss' 
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(10) 
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is the most probable velocity in a gas of particles of mass jj^ 
and temperature 



-* ss' — 
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The quantity {crv}ss' is called the momentum transfer rate coeffi- 
cient. A related quantity is the collision frequency between a test 
particle s and the incident particles s' 



Vss' = ns'{o-v}s 



(13) 



The momentum transfer rate coefficient defined by eq. (|8j is a 
function of the temperature Tss' and the drift velocity |Uj' - Uj|, 
with the parameters x and ^ playing the roles of Mach number 



for the collision and the drift velocity, respectively. Owing to the 
dependence of {crv)ss' on ^, the friction force eq. (|6]l is in general 
a nonlinear function of the drift velocity. 

For a cross section varying as a power-law of the relative ve- 
locity, <t(Vss') = o-qx^", the integral in eq. dS) can be expressed 
in terms of the confluent hypergeometric function M(a,b,c) 
(Abramowitz & Stegun 1965), 

<crv)„ - croass'Gni^), (14) 

where 

For vanishing drift velocities G„(^) reduces to a constant factor, 
G„(0) = 4r(3 - n/2)/(3 y/n), whereas in the opposite limit of 
large drift velocities G„(^) ~ ^'"" . The behaviour of the function 
G„(^) for n - 0,1,2,3 and 4 is shown in Fig.[T] 

Three special cases of eq. ( fT4l ) are relevant for collisions be- 
tween ISM particles: 

(a) « = 0, with 



""^^'^'-T^'^wy'^i'^r^^''''''' 



(16) 



(Baines et al. 1965). A cross section ctq - n(rs -H rs')^, inde- 
pendent on the relative velocity describes the collisions between 
solid particles, represented as "hard spheres" of radii r^ and r,/ . 

(b) n - \, with Gi(^) = 1. A cross section varying as 
the inverse of the relative velocity arises for an induced dipole 
attraction, where the interaction potential is proportional to 
the inverse fourth power of the distance ("polarization poten- 
tial"). Particles obeying this particular interaction law are called 
"Maxwell molecules"(Maxwell 1860a,b), and the corresponding 
collisional rate coefficient "Langevin rate" (see Appendix A). In 
this polarization approximation, the momentum transfer rate is 
independent of both temperature and drift velocity. 

(c) n - 4, corresponding to a (screened) Coulomb potential, 
with 






2^ 



(17) 



This result was first derived by Chandrasekhar (1943) in the con- 
text of stellar dynamics and, independently, by Dreicer (1959) 
and Sivukhin (1966) for a fully-ionized gao 

In many astrophysical applications, the relative velocity of 
the interacting particles is of the order of the sound speed in the 
gas (if the temperatures of the two species are not too differ- 
ent), but the drift velocity lu,- - u,! is usually much lower than 
the sound speed. In such circumstances it is appropriate to ex- 
pand the integrand in eq. ^ for small values of the quantity ^x, 
obtaining, to the lowest order in the expansion, a value of the 
momentum transfer rate coefficient independent of the drift ve- 
locity. 
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' Draine & Salpeter (1979) give the following accurate approxima- 
tions f or « = and n = 4: 
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Fig. 1. The normalized function G„(^)/G„(0) defined by eq. ( fTsT i 
as function of the normalized drift velocity ^ = lu, - u.s'|/fl.5j' 
with « = 0, 1, 2, 3, and 4 (solid curves). The dashed curves 
show the approximations given by Draine & Salpeter (1979) for 
the cases n - and n - A. 



If the momentum transfer cross section is given in terms of the 
collision energy in the center-of-mass system £„„ = fj.ss'v\^,/2, 
the rate coefficient in the small drift limit can be written as 



<0-v)s 
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(19) 



where z = E^m/kTss', a frequently-adopted expression (see e.g., 
Mitchner & Kruger 1973).B 



3. Collisions with H2 

Collisions between H2 molecules and charged particles deter- 
mine the rate of diffusion of the interstellar magnetic field 
through the dominantly neutral gas of cold molecular clouds (see 
Paper I). In applications to molecular clouds (e.g., Nakano 1984, 
Mouschovias 1996), the collision rate coefficient {crv}s^H^ has 
usually been estimated with the polarization approximation 
(see Appendix) for collisions with molecular ions, and, also 
sometimes, for collisions with electrons (see discussion in 
Mouschovias 1996). For collisions between H2 and dust grains, 
the hard sphere model (Sect. 2) has generally been assumed. We 
review below the validity of these assumptions and compute ac- 
curate values for the collision rate coefficients using available 
momentum transfer cross sections. 



^ The momentum transfer cross section measured in laboratory ex- 
periments is usually given as function of the energy in the laboratory 
frame (£iab)- Before using eq. ( |19l >, iiiab must be converted into the en- 
ergy in the center-of-mass frame according to the formula 
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where s and s' indicate the target and incident particles, respectively. 



Fig. 2. The momentum transfer cross section for collisions 
HCO^-H2 as a function of £„„ computed by Flower (2000) (dots 
and solid curve). The dashed curve shows the Langevin cross 
section. 



3.1. HC0^-H2 

Flower (2000) has calculated quantum-mechanically the cross 
section between HCO^, a dominant ion in typical molec- 
ular cloud conditions, and H2 molecules in their rotational 
ground states. Figure |2] shows the cross section computed by 
Flower (2000) compared with the Langevin cross section. As 
noticed by Flower (2000), the Langevin value gives a good ap- 
proximation to the quantal results. The rate coefficient is, there- 
fore, expected to depend very weakly on temperature and drift 
velocity, as shown in Fig. (Oq Previous estimates of the molec- 
ular ion-H2 rate coefficient based on the Langevin formula (see 
Fig. [3]) have ignored the weak dependence from temperature and 
drift velocity. This neglect is, however, of little consequence for 
models of magnetically-controlled cloud collapse as long as the 
temperature of the infalling gas is in the range 10-20 K and the 
ion-neutral drift velocity is a small fraction of the sound speed, 
as shown e.g., by Mouschovias & Ciolek (1994). 



3.2. H+ - H2 

Figure |4] shows the momentum transfer cross section for colli- 
sions of H3 and H2, recommended by Phelps (1990), converted 
to energies in the center-of-mass frame. At low energies (E^m ^ 
1 eV), the cross section is obtained from the results of the mo- 
bility experiments of Ellis et al. (1976). At higher energies, the 
behavior of the cross section is extrapolated from mobility data 
at room temperature. As shown in the figure, the cross section is 
very close to the Langevin value up to E^m ~ 3 eV, but declines 
steeply (as E^^'') above Ecm ~ 10 eV. Since the cross section tab- 
ulated by Phelphs (1990) does not extend to ion energies below 



' Our rate coefficient for zero drift velocity differs by ~ 20% from 
the result by Flower (2000) who adopts an averaging over the particles 
kinetic energy slightly different from our eq. ( 1191 ). 
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Fig. 3. The rate coefficient for HCO^-H2 collisions as a func- 
tion of the temperature T for Vd - (thick solid curve); and 
Vd - 1 km s"'; Vd - 2 km s"' (thin solid curves, bottom to top), 
compared with the Langevin rate (dashed line). 



Fig. 4. Momentum transfer cross section for H^ -H2 collisions as 
a function of the energy in the center-of-mass system according 
to Phelps (1990) (solid curve). The dotted line shows our ex- 
trapolation at low energies, whereas the dashed line shows the 
Langevin cross section. 



0.1 eV, we have extrapolated Phelps's value to lower energies 
with the asymptotic formula cr^i - 1.62 x 10"'-''(£'cm/eV)""-''°^. 
Figure |5] shows the corresponding rate coefficient obtained inte- 
grating the cross section, according to Eq. ([8j as a function of the 
temperature and for various values of the drift velocity. As in the 
case of HCO^, given the Langevin behavior of the cross section 
at low energies, the dependence of the collisional rate from tem- 
perature and drift velocity is very weak for temperatures below 
~ lO'* K. 



3.3. H+-H2 

Figure |6] shows the momentum transfer cross section for H^- 
H2 collisions computed by Krstic & Schultz (1999) in the en- 
ergy range 0.1 eV < E^n < 100 eV (fully-quantal calcu- 
lation), by Bachmann & Reiter (1995) in the energy range 
0.05 eV < /icm < 10^ eV (classical calculation) and the cross 
section recommended by Phelps (1990) in the energy range 
0.067 eV < ^cm < 6.67 keV (obtained interpolating results 
from ion mobility experiments for /Scm < 2.3 eV and theoretical 
cross sections above 330 eV). As shown in Fig. |6] the quantal 
results depart significantly from the Langevin value at high en- 
ergies (£■„,! > 4 eV), and differ substancially from the results 
extrapolated from ion mobility experiments in the same energy 
range. Here we adopt the cross section computed by Krstic & 
Schultz (1999), extrapolated to energies below 0.1 eV with the 
asymptotic formula cr^i - 1.05 x 10"'^(£'cm/eV)"''-^^'. We show 
in Fig. I2] the corresponding rate coefficient as a function of the 
temperature and the drift velocity. Again, the rate coefficient de- 
pends very weakly on these two quantities for temperatures be- 
low ~ 10^ K. 



3.4. e-H2 

The elastic scattering of electrons by H2 has been the subject 
of much theoretical and experimental work of the past decades. 
In general, it is well established that the momentum transfer 
cross section for electron-molecule scattering deviates signifi- 
cantly from the classical Langevin rate at low energies (below 
~ 1 eV), owing to the effect of "electron exchange", i.e., the 
exchange of the incoming electron with one orbital electron in 
the neutral species (see e.g., Massey & Ridley 1956, Morrison 
& Lane 1975). The existence of this effect has been confirmed 
experimentally by Ferch et al. (1980) in the energy range 0.02- 
2 eV. The net result is a reduction of the momentum transfer 
cross section, which compensates the polarization (Langevin) 
contribution in such a way that the resulting cross section is 
roughly constant at low-collision energies. 

We have assembled a compilation of the available measure- 
ments of the e-H2 momentum transfer cross section, for collision 
energies ranging from 10"^ eV up to 200 eV. A detailed sum- 
mary of the most recent experimental results is given by Brunger 
& Buckman (2002). The results are shown in Fig. [8] compared 
with the theoretical calculations of Henry & Lane (1969) and 
the Langevin value. Clearly, the semi-classical Langevin for- 
mula provides a poor approximation to the actual cross section, 
especially at energies below ~ 1 eV, where the effects of "elec- 
tron exchange" are dominant. The agreement between labora- 
tory measurements and theoretical values of the scattering cross 
sections appear to be satisfactory, with the possible exception 
of the region of very low electron energies, where the possible 
existence of a Ramsauer-Townsend minimum in the momentum 
transfer cross section is not completely excluded (Ramanan & 
Freeman 1991). The momentum transfer collision rate calcu- 
lated with the cross section indicated by the solid line in Fig. |8] 
is shown in Fig.|9] At temperatures typical of interstellar clouds. 
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Fig. 5. The momentum transfer rate coefficient for H3 -H2 col- 
lisions as function of the temperature T computed with the 
cross section shown in Fig. |4] for vj = (thick solid curve); 
Vd - I km s"' (dotted curve); v^ = 5 km s"' (short-dash dot- 
ted curve); and Vd - 10 km s"' (long-dash dotted curve). The 
dashed line shows the Langevin rate. 



Fig. 7. The collision rate coefficient for H^-H2 collisions as 
a function of the temperature, computed with the momentum 
transfer cross section of Krstic&Schultz( 1999) for Vrf = Q(thick 
solid curve); and v^ = 2 km s"'; Vd - 4- km s"'; Vd = 10 km s"' 
(thin solid curves, bottom to top). The dashed line shows the 
Langevin rate. 
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Fig. 6. The momentum transfer cross section for H^-Hi colli- 
sions as a function of E^m'- fuUy-quantal calculation by Krstic & 
Schultz (1999) (solid curve), with our extrapolation to low ener- 
gies (dotted line); classical calculation by Bachmann & Reiter 
(1995) (long-dashed curve); semi-empirical recommendations 
by Phelps (1990) (dot-dashed curve); Langevin cross section 
(dashed curve). 



T X 10 K, the Langevin formula overestimates the actual value 
of the momentum transfer rate by about two orders of magni- 
tude. 



4. Collisions with H 

4.1. C+ -H 

The momentum transfer cross section for collisions of C^ 
ions with H atoms was computed by Flower & Pineau-des- 
Forets (1995) at several values of the collision energy, ranging 
from Eon ~ 5 X lO"'' eV to £„„ a; 5 eV with the adiabatic po- 
tential of Green et al. (1972) (see Fig.fTOji. Figure [TT] shows the 
corresponding rate coefficient computed by integrating numer- 
ically our interpolation of the results of Flower & Pineau-des- 
Forets (1995). 

4.2. H+ -H 

Being the simplest ion-neutral collision process, the scattering 
of H^ by H atoms has been the subject of a large number of 
theoretical investigations, from the semiclassical calculations of 
Dalgarno & Yadav (1953) to the accurate quantum mechanical 
calculations of Krstic & Schultz (1999) and Glassgold, Krstic 
& Schultz (2005, hereafter GKS). The latter paper, reporting a 
determination of the momentum transfer cross section over a 
range from 10"'*' eV to 10^ eV (partially shown in Fig. fT2l i. 
represents the definitive reference on this elastic process. The 
only experimental determination of the H-H^ momentum trans- 
fer cross section was obtained by Brennan & Morrow (1971) at 
Eon ~ 5 eV by observing the velocity and attenuation of com- 
pressional Alfven waves propagating in a partially-ionized hy- 
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Fig. 8. The momentum transfer cross section for e-H2 colli- 
sions as a function of the electron kinetic energy. Experimental 
values: England et al. (1988) (empty triangles), Schmidt et 
al. (1994) (empty squares), Shyn& Sharp (1981) (empty circles), 
Nishimura, Danjo & Sugahara (1985) (filled triangles), Khakoo 
& Trajmai- (1986) (filled squares), Brunger et al. (1990,1991) 
(filled circles). The dashed and dotted curves show the Langevin 
cross section and the quantum-mechanical theoretical results of 
Henry & Lane (1969), respectively. The solid curve shows the 
cross section adopted in this work. 



drogen plasma. The measured value, 0-^^,^ - 6.0^j5 x 10"'^ cm^, 
is in good agreement with the quantum mechanical results (see 
Fig. [T2I 1. Previous estimates of the momentum transfer rate co- 
efficient (Geiss & Biirgi 1986) were based on measurements 
and calculations of the charge exchange cross section a-^e for 
the reaction H -H H^ — » H^ H- H complemented at low energies 
with the polarization cross section computed with the Langevin 
formula eq. jAAj . Since the charge transfer process proceeds 
with little momentum transfer between the interacting particles, 
eq. (HJi with ® - n gives 0-^1 ~ ^a-^.^ (Dalgarno 1958, Banks 
& Holzer 1968). As shown by Fig. [121 ^ combination of the po- 
larization cross section and twice the value of the charge ex- 
change cross section roughly reproduces the accurate results of 
GKS. However, the momentum transfer rate estimated by Geiss 
& Biirgi (1986) with this approximation for temperatures be- 
tween 10^ K and 2 x lO'* K is about 50% higher than the rate 
obtained by a numerical integration of the momentum transfer 
cross section of GKS, shown in Fig. [T3] for different values of 
the drift velocity. 



4.3. e-H 

The momentum transfer cross section for e-H collisions has 
been computed by Dalgarno, Yan, & Liu (1999) with the quan- 
tal formula of Dalgarno & Griffing (1958) and the phase shifts 
of Rudge (1975) and Das & Rudge (1976) at £ < 10 eV, and 
by van Wyngaarden & Walters (1986) in the energy range 100- 
300 eV. Laboratory measurements of the e-H momentum trans- 



Fig. 9. The momentum transfer rate coefficient for e-H2 colli- 
sions computed with the cross section shown in Figl8]as a func- 
tion of the temperature T for v^ = (thick solid curve); and 
Vd - 50 km s"'; vj = 100 km s"' (thin solid curves, top to 
bottom), compared with the values computed with Eq. (14) of 
Draine (1983) for the same values of the drift velocity (dotted 
curve).The dashed line shows the Langevin rate. 




Fig. 10. The momentum transfer cross section for collisions C^- 
H computed by Flower & Pineau-des-Forets (1995) as a func- 
tion of the collision energy in the center-of-mass frame (original 
data, triangles; our interpolation, solid curve). The dashed curve 
shows the Langevin cross section. 
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Fig. 11. The collision rate coefficient for collisions C^-H as a 
function of the temperature for vj - (thick solid curve); and 
Vd = 5; v^ = 10; and v^ = 15 km s"' (thin solid curves, bottom 
to top), compared with the Langevin rate (dashed line). 
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Fig. 12. The momentum transfer cross section for collisions H^- 
H computed by GKS as a function of the collision energy in the 
center-of-mass frame (solid curve), compared to the Langevin 
cross section (dashed curve). The empty circle is the measure- 
ment by Brennan & Morrow (1971) at E^m ~ 5 eV. The filled tri- 
angles are the values of the charge exchange cross section for the 
same interacting particles, measured by Newman et al. (1982), 
multiplied by 2. 



Fig. 13. The collision rate coefficient for collisions H^-H as a 
function of the temperature for vv = (thick, solid curve); vj = 
10; Vd - 50 km s"'; and vj - 100 km s"' (thin, solid curves, 
bottom to top); compared with the fitting formula of GKS for 
the same values of the drift velocity (dotted curves), the fitting 
formula of Geiss & Biirgi (1986) for v^ = (long-dashed curve) 
and the Langevin value (short-dashed curves) . 



fer cross sections have been performed by Williams (1975a,b), 
Callaway & Williams (1975), Shyn & Cho (1989), and Shyn & 
Grafe (1992) (for a detailed review of the experimental meth- 
ods and results, see Bederson & Kieffer 1971 and Trajmar & 
Kanik 1995). Theoretical and experimental results are shown in 
Fig. [14] The dependence of the cross section on collision en- 
ergy is clearly non-Langevin: at low energies the cross section 
is approximately constant, cr k 4 x 10"'^ cm"^, whereas at high 
energies, it decreases with energy as E^^^. As a result, the mo- 
mentum transfer rate, shown in Fig. [15] is lower by about one 
order of magnitude than the Langevin value at T x 10 K, and 
larger by ~ 40% at T x 10'* K, with a weak dependence on the 
relative drift velocity. 



5. Collisions with He 

With the exception of collisions with H^ or electrons, to 
our knowledge, there are no experimental or theoretical data 
available for elastic collisions between He atoms and charged 
species. An estimate based on the Langevin approximation (see 
Appendix), gives 



<0"V),-He = 



(m,- H- miie)m„piie 



(nij + m„)miiep„ 



1/2 



{o-v}i, 



(20) 



where n -H or H2, and p„ is the polarizability of species n (see 
Appendix). According to this equation, the momentum transfer 
rate for collision of ions with He is, in general, a factor 0.4-0.5 of 
the corresponding rate for collisions with H2 and a factor 0.3-0.4 
of the corresponding rate for collisions with H. 
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Fig. 14. Momentum transfer cross section for e-H collisions. 
Experimental values: Shyn & Cho (1989) (filled circles); 
Williams (1975a,b) (filled triangles); Shyn & Grafe (1992) 
(filled squares); Callaway & Williams (1975) (empty circles). 
The solid curve is the theoretical calculation of Dalgarno et 
al. (1999). 



Collisions with He introduce a small correction to the ex- 
pression of the friction coefficient a,,, defined by eq. O. For a 
neutral component made of H2 and He, eq. ( l20b gives 



a in = ff/.H, + O'l.He - C/a/.Hj, 



with 



1 + 



(nii + mH,)miiePHe 



(nii +mHe)wH,;?H2 



1/2 



/ "He 



(21) 



(22) 



For example, for a cosmic He abundance of 0.1, the He 
correction factors based on the Langevin approximation are 
ch+ - 1.12, ch+ = 1.13, and Chco+ - 1-15 (see also 
Mouschovias 1996). Similarly, for collisions with atomic hydro- 
gen, the He correction for the H-H^ rate coefficient is ch+ - 
1.08. 



5.1. H^ -He 

Figure [16] shows the momentum transfer cross section for H^- 
He collisions computed by Krstic & Schultz (1999) with a semi- 
classical treatment in the energy range 0.1 eV < E^^ < 100 eV. 
The calculation has been recently extended up to 10"^ eV by 
Krstic & Schultz (2006). The theoretical results are in good 
agreement with the Langevin value below Ecm ~ 1 eV, but at 
larger energies the Langevin formula overestimates the theoret- 
ical results (by one order of magnitude at Zscm ~ 100 eV). To 
our knowledge, no experimental results are available for H^- 
He collisions. Here we adopt the cross section computed by 
Krstic & Schultz (1999), extrapolated to energies below 0.1 eV 
with the Langevin value. The resulting collisional rate coeffi- 
cient is shown in Fig. [T7] as a function of the temperature and 



Fig. 15. Collision rate coefficient for e-H collisions as a function 
of temperature, for v,/ = km s"' (thick curve); vj - 50 km s"'; 
and Vrf - 100 km s"' (thin curves, bottom to top). The dashed 
line shows the Langevin rate. 



the drift velocity. As for H^-H2 collisions, the rate coefficient 
depends very weakly on these two quantities for temperatures 
below ~ lO"* K and drift speeds below ~ 10 km s"' . The rate co- 
efficient (crv)H+,He IS about 0.9 (crv)H+,H2J whereas eq. (|20] i gives 
a factor 0.5. 



5.2. e-He 

Accurate values of the momentum transfer cross sections for 
e-He collisions have been obtained from mobility experi- 
ments (Crompton, Elford & Jory 1967; Crompton, Elford 
& Robertson 1970; Milloy & Crompton 1977; Ramanan & 
Freeman 1990). The agreement between different experimental 
evaluations is excellent, at the level of 1-2%. As in the case of 
e-H2 collisions discussed in Sect. [3] the actual cross section de- 
viates significantly from the classical Langevin value, owing to 
the quantum exchange of the incoming electron with one orbital 
electron of He. This is especially evident at low-collision en- 
ergies, where the momentum transfer cross section is approxi- 
mately constant. The rate coefficient (crv)e_He is a factor ~ 0.6 of 

<0"v)e,H,- 



6. Collisions between dust grains and neutral 
particles 

For collisions between a spherical grain g of charge ZgC and ra- 
dius Tg with a neutral particle n of polarizability a„, the geomet- 
rical cross section nr^ is larger that the polarization cross section 



(eq. IA.4I ) when 



rg > 24|Zj 



,1/4 



aV UJ 



-1/4 



(23) 



Pinto & Galli: Momentum transfer rate coefficients 



10-14 



10- 



10- 



10- 



0.01 




100 




Fig. 16. The momentum transfer cross section for collisions 
H^-He, according to the fully-quantal calculations of Krstic & 
Schultz (1999), as a function of the collision energy in the center- 
of-mass frame (solid curve), compared to the Langevin cross 
section (dashed line). 



Fig. 17. The collisional rate coefficient for H^-He collisions as 
a function of the temperature for v,/ = (thick, solid curve); 
Vd = 10 km s"'; v^ = 20 km s"'; and v^ = 30 km s"' (thin, 
solid curves, top to bottom) compared with the fitting formula 
of Geiss & Biirgi (1986) for v^ = (long-dashed curve) and the 
Langevin value (short- dashed curves). 



a condition generally satisfied by interstellar grains. The aver- 
age momentum transfer cross section for collisions of spherical 
grains with neutral atoms or molecules is then 



<crv), 



gn 






(24) 



where Go(^) is given by eq. (fTSI l and 5 is a factor of order unity 
(the so-called Epstein coefficient) equal to unity if the neutrals 
impinging on the grain undergo specular reflections. For very 
subsonic drift velocity, this expression reduces to 



<crv)g„ * -Ttr^d 



UT„ 



1/2 



(25) 



a result derived by Epstein (1924). Experiments with micron- 
size melamine-formaldehyde spheres show that 6 ~ 1.3 (Liu et 
al. 2003). 

Figure |20] shows the grain-neutral momentum transfer rate 
as a function of the relative drift velocity according to eq. jTM 
with (5 = 1 compared with the approximations given by 
Draine & Salpeter (1979), Nakano (1984), and Mouschovias & 
Ciolek (1999). In this figure, the grain-neutral momentum trans- 
fer rate is normalized to the expression given by Mouschovias 
& Ciolek (1999) in the low-drift limit, and the drift velocity is 
normalized to the mean thermal speed in the neutrals. 



7. Collisions between charged particles 

Particles of charge Z^/g can exchange momentum with parti- 
cles of charge Zjg and density n, via long-range Coulomb in- 
teractions. In the standard Coulomb scattering theory (see e.g.. 




100 



Fig. 18. The momentum transfer cross section for e-He colli- 
sions as a function of the electron kinetic energy. Experimental 
values: Crompton et al. (1970) (filled triangles), Milloy & 
Crompton (1977) (filled squares). The dotted curve shows the 
experimental results of Ramanan & Freeman (1990), with un- 
certainty ~ 3%. The dashed curve shows the Langevin cross 
section, whereas the solid curve is the cross section adopted in 
this work. 
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Fig. 19. The momentum transfer rate coefficient for e-He col- 
lisions computed with the cross section shown in Fig[T8] as a 
function of the temperature T for v^ = (thick, solid curve); and 
Vd = 50 km s"'; Vd = 100 km s"' (thin, solid curves, bottom to 
top). The dashed Une shows the Langevin rate. 



Fig. 20. Momentum transfer rate coefficient for collisions grain- 
neutral as a function of the relative drift velocity, from eq. ( fTSl ) 
with n = (solid curve), compared with the expressions from 
Draine & Salpeter (1979) (dotted curwe.); Nakano (1984) (dashed 
curve); and Mouschovias & Ciolek (dot-dashed curve). 



Chapman & Cowling 1953), the momentum transfer cross sec 
tion is given by 



CTmt - 27rrc(v„') In 



r2 

max 



'"c^^s^ 



) + r 



,2 
min 



where 



rc(vss') 



ZjcZjc' € 



(26) 



(27) 



is the Coulomb radius (the distance at which the electrostatic en- 
ergy is of the order of the kinetic energy) at the relative velocity 
Vss', and rmiii, r^ax are the minimum and maximum impact pa- 
rameters of the collision. Several different approaches have been 
taken to calculate r^in and r^^^. For the latter, the standard choice 
is the total Debye length in the plasma Ad (Cohen et al. 1950), 
given by 

D s D,s s 

When averaging over the relative velocity of interacting par- 
ticles (eq. Is), the slowly varying logarithmic factor can be con- 
sidered constant, and the term rc(vss') in the logarithm's argu- 
ment can be replaced by 



re = rc({v,,>}) = 



3kT,. 



(29) 



where (v^^,) = 3kTss'/tJ^ss' is the average of the square of the 
relative velocity (Spitzer 1956). The averaging of the remaining 
dependence of the cross section on the inverse fourth power of 
the relative velocity then gives, according to eq. ( fTSl ). 



(crv)jj' = 27:ass'rr(ass')in 



W + 4 



re 



G4(0, 



(30) 



where ass' is given by eq. ( fTTT i and G4(^) by eq. ( fTTI ). If the par- 
ticles are assumed to be pointlike (rmi„ = 0), and, in addition, 
^D » /t (a condition generally fulfilled in the ISM), we have 



In 



1+1^ 
re 



21n|^i. 

re. 



(31) 



In this case eq. ( l30b for ^ = gives the standard expression of 
the Coulomb momentum transfer rate (Spitzer 1956) 

len'/^Z^Zle' 
<o"v),y = — — --rlnAi.5' 

3i^lih2kTss')''^ 

.8.48xlO-2zX^,(^)"'(|pinA, (32) 



where 



InAtt' = ln| ^1 = In 
re 



3kTss_(4ne^y^ns] ^'^ 



9.42 - In < 



s s' \ 

z2(«,/cm-3) 



Ts 



^(Tss'l^) 

is the Coulomb logarithm. 



E 



(r.,/K) 



(33) 
(34) 



7.1. Collisions of ions with charged dust grains 

In astrophysical applications, eq. ( |32] ) has often been adopted to 
compute the rate of momentum transfer between grain-ion and 
grain-electron collisions (Draine & Salpeter 1979). However, 
a more accurate calculation requires a modification of the 
Coulomb logarithm to account for the finite size of the grains. 
Taking rmin = rg in Eq. ( [30l l. and assuming as before Au » re. 
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- = rkV(Zje)= 



Fig. 21. Argument of the Coulomb logarithm A^, for grain-ions 
interactions normalized to Af,- for pointlike particles, and ratio 
fglTc as a function of the reduced temperature of the ions r, - 
rgkTj/iZje)^. The solid and dashed curves are for "light ions" 
with effective atomic weight 1 and "heavy ions" with effective 
atomic weight 25. 



Table 1. Fitting formulae for momentum transfer coefficients as 
a function of the gas temperature T (in K) and 6 = log(r/K) for 



species s, s' (crv)ss' 
(10-" , 



-') 



HCO+,H2 

H3,H2 

H+,H2 
e.H, 



r"-(1.476 - 1.4096/ + 0.5550- - 0.077561') 
2.693 - 1.2380 + 0.6646*2 - 0.0896*^^ 
1.003 + 0.0500 + 0.13602 - 0.0140^ 
r'/2(o.535 + 0.2030 - 0.16302 + 0.0500^) 



C+,H 
H+,H 



1.983 + 0.4250 - 0.43102 + 0.1140^ 
0.6497-0.375 (Qj^s) 

r'/2(2.841 + 0.0930 + 0.24502 - 0.0890^^) 



H+,He 
e,He 



1.424 + 7.438 x lO^T - 6.734 x 10-^2 
0.428r'/2 



Table 2. Fitting formulae for momentum transfer coefficients as 



function of v'rms (in km s ). 



species s, s' (crv),,/ 

(cm^^ s"') 



HCO+,H2 

H3,H2 

H+,H2 



2.40 X 10-%i'„^, 
2.00 X 10-'v°„'5 

3.89 X lO-'v.-r' 
3.16 X 10-"v', 



(kms-') 



rms 



0.2 < VVms < 5 
1 < Vn,s < 10 
1 < V„,, < 10 

20 < v,.„,s < 200 



C+,H 
H+,H 
e,H 



2.13x10-^,75 
2.50xl0-"'vv',^,exp(-v„,,/460) 



2 < Vn,, < 20 

v„„, > 1 (GKS) 
20 < v,,^s < 600 



H+,He 
e,He 



1.48 X 10- 
7.08 X 10- 



0.1 <vv„,, < 10 
20 < v,,„s < 500 



we obtain the modified Coulomb logarithm for grain-ion (or 
grain-electron) collisions 



lnA„ 



In 



(7t' + rp'/2 



(35) 



which reduces to ln{A£,/7c) for r^ - 0, and approaches \n(A£,/rg) 
for Tg » Tq. The latter value of the Coulomb logarithm has been 
used, e.g., by Benkadda et al. (1996) in their study of nonlinear 
instabilities in dusty plasmas. Following Draine & Sutin (1987), 
it is easy to compute the ratio rg/Tc for a variety of astrophysical 
conditions, since this quantity depends almost exclusively on the 
reduced temperature T/ = rgkTiliZie)^, 



re 



3Z,T,- 

ZgiTiY 



(36) 



Using the analytical approximations for the value of the aver- 
age grain charge Z^e as a function of t, computed by Draine 
& Sutin (1987), we obtain the modified Coulomb logarithm for 
grain-ion collisions given in Fig.|2T|for "light" and "heavy ions" 
(effective atomic weight 1 and 25, respectively). We see from 
this figure that the pointlike approximation for the Coulomb log- 
arithm remains valid for t, less than ~ 0. 1 . For larger values of r, 
the value of the Coulomb logarithm is smaller than the pointlike 
value, but the effect is a reduction of the constant in Eq. ( l34b . 
from 9.42 to 8.97 or 9. 18 at worst, for the "light ion" and "heavy 
ion" cases, respectively. 

For a more general calculation of the Coulomb logarithm for 
grain-ion and grain-electron where r^i^ depends on the finite size 
of the grain and the relative velocity of the collision, and, in ad- 
dition, the approximation Ao '» Tc ^& removed, see Khrapak & 
Morfill (2004). 



8. Analytical approximations 

In Table [T] we list analytical fitting formulae for the momentum 
transfer rate coefficients computed numerically in Sects. 3-5. 
For zero drift velocity, we have approximated the rates (or their 
logarithm) with third-order polynomials of the logarithm of the 
temperature. The accuracy of these fitting formulae is ~ 10%. 
The dependence on the drift velocity has been approximated, 
following GKS and Draine (1980) by power-laws (or modified 
power-laws) fits of the numerical results in terms of the rms ve- 
locity, defined as 



^'d + — - 



TTf^s 



1/2 



(37) 



These fitting formulae are shown in Table |2] Within the range of 
Vims indicated, the accuracy of the power-law approximations is 
better than ~ 20%. 



9. Conclusions 

We have derived momentum transfer coefficients for col- 
lisions between ions, electrons, charged dust grains and 
atomic/molecular hydrogen from available experimental data 
and theoretical calculations, within the classic approach devel- 
oped by Boltzmann (1896) and Langevin (1905). The numerical 
results have been approximated with simple analytical functions 
of the temperature and the drift velocity between the colliding 
species. The main conclusions of our study are the following: 

1. For collisions between molecular ions (HCO^, Hp and H2, 
the often used polarization approximation valid for an in- 
duced electric dipole attraction is in satisfactory agreement 
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(within 20-50% of the numerical results). At low tempera- 
tures (10-50 K) and for low values of the drift velocities, the 
polarization approximation overestimates the actual rate for 
collisions between H^ and H2 by a factor of ~ 3. In heavily- 
depleted molecular cloud cores, where H^ may be the dom- 
inant ion (see, e.g., Walmsley et al. 2004), the Langevin ap- 
proximation should not be used to compute ambipolar diffu- 
sion time scales (see Paper I). 

2. For collisions between electrons and H2, which contribute 
marginally to the resistivity of the interstellar gas (see 
Paper I), the polarization approximation fails by orders of 
magnitude because of the non-dipolar nature of the interac- 
tion, especially at low energies. 

3. In the range 10K< T < 10^ K, rate coefficients for colli- 
sions between H^ and He atoms are a factor 0.7-1 and 0.1- 
0.4 of the corresponding rate for collisions with H2 and H, 
respectively. In the same temperature range, the rate coeffi- 
cient for e-He collisions is a factor 0.4-0.7 of the rate for 
e-H2 collisions and a factor 0. 1-0.3 of the rate for e-H col- 
lisions . 

4. Grain-neutral and grain-ion collisions are well represented 
by hard-spheres and Coulomb interactions, respectively. In 
the latter case, slightly different expressions may be adopted 
for the Coulomb logarithm, with no significant consequences 
for applications to ISM grains. 
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supported by the Marie Curie Research Training networks "Constellation". 



Appendix A: Polarization approximation 

The cross section for the interaction between an ion i of charge Z,e and a neutral 
molecule (or atom) n is determined by the attractive polarization potential 



V,n(.r) ■- 



2r* 



(A.l) 



where r is the distance between the centers of the ion and the molecule (much 
larger than the sizes of the ion and molecule) and a„ is the polarizability 
of the molecule. The relative cross section can be calculated from the clas- 
sical trajectories allowed by this potential (see, e.g., Dalgamo, McDowell & 
Williams 1958a,b; Mitchner & Kruger 1973), and is 



o-(v,„) = 2.210;t 






1/2 



(A.2) 



In this polarization approximation, the rate coefficient ("Langevin rate") is given 
by 



<o-v>,„ = 2.2l0n 



Pn7.se 



,2\l/2 



(A.3) 



which is independent of ffie gas temperature and the relative drift of ffie interact- 
ing species. Osterbrock (1961) corrected the numerical coefficient 2.21 in 2.41 to 
account for ffie repulsive nature of the interaction potential at small r, assuming 
that for small impact parameters the scattering is on average isotropic. 

The Langevin cross section for collisions of ions with neutrals, taking into 
account Osterbrok's (1961) correction, therefore, results in 



^,„(£„,) = 2.()3xlO-'5z'/2h?^ {^] cm2. 



1/2, p ,-1/2 



(A.4) 



and the corresponding rate is 



1/2 , v-l/2 



<o-v),-, 



^2.S1X10-Z-(||) (^) 



(A.5) 



Values of ffie polarizability are pn = 0.667 A^, pn, = 0.804 A^, and pHe 
0.207 (Osterbrock 1961). 
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